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Abstract. We discuss the photoionization of the inter- 
galactic medium by early formed quasars and propose 
a mechanism which could generate measurable temper- 
ature fluctuations 6t = AT/T of the cosmic microwave 
background (CMB). Early quasars produce individual ion- 
ized regions around themselves. We evaluate both ther- 
mal and kinetic (first-order Doppler) Sunyaev-Zeldovich 
(SZ) effects associated with those ionized regions. Whereas 
the former is negligible, the latter induces the generation 
of detectable secondary small and medium scale CMB 
anisotropies. We find that this effect could produce mea- 
surable individual sources with a St — 10~^ to ~ 10~* on 
scales up to 1°, which introduce a non gaussian signature 
in the statistics of the primordial anisotropies. The effect 
of the early formed quasars is compatible with the sources 
detected in Cheng et al. (1994) MSAM CMB experiment. 

Key words: cosmology: cosmic microwave background - 
galaxies: quasars 



1. Introduction 

In most cosmological models, some reionization of the in- 
tergalactic medium (IGM) must have occur between re- 
combination and the present time where most of the gas 
is known to be already reionized (Gunn & Peterson 1965). 

Several mechanisms were used to explain the ionization 
of the IGM. The possibility of photoionization has been 
extensively explored (Shapiro & Giroux 1987; Donahue & 
ShuU 1987; Miralda-Escude & Ostriker 1990; 1992). The 
possible sources of photons like young galaxies, primordial 
stars and decaying particles do not seem to be entirely 
satisfactory. 

Apart from photoionization, the heating of the IGM 
may occur in some scenarios of galaxy formation involv- 
ing explosions (Ikeuchi 1981; Vishniac & Ostriker 1985). 
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Another possibility, which arises in hierarchical models, is 
ionization by high mass stars which form in subgalactic 
size clumps at high redshifts (Couchman & Rees 1986). 

A detailed discussion of the various processes of reion- 
ization would be beyond the scope of this paper. We will 
focus our interest on the reionization by early quasars and 
how it affects the cosmic microwave background (CMB), 
knowing that Meiksin & Madau recently (1993) have 
shown that the quasars might indeed produce enough ion- 
izing photons to completely ionize the IGM. 

The main effect of the reionization is to wash out some 
of the temperature fluctuations of the CMB on small scales 
which undergo Thomson scattering on the reionized IGM 
(see for example Bond & Efstathiou 1984; Vittorio & Silk 
1984; Tegmark et al. 1994). 

In the assumption of a partial reionization by quasars, 
we present a new effect which, on the contrary, generates 
secondary small scale distortions on the CMB. 

At times preceding the full reionization, quasars ionize 
the gas around themselves and they become surrounded 
by ionized bubbles. Those bubbles are moving together 
with the quasars in the large scale gravitational potential 
wells. Thus they have got a bulk peculiar velocity v, with 
respect to the local standard of rest in which the CMB is 
isotropic. 

An effect equivalent to the Sunyaev-Zeldovich (SZ) ki- 
netic effect for rich clusters creates relative temperature 
fluctuations which are proportional to the product of the 
Thomson optical depth t by the ratio of the quasar radial 
velocity to the velocity of light (u, /c) - between a line of 
sight crossing the ionized bubble and a nearby one which 
does not. Contrary to X-ray clusters , the corresponding 
thermal SZ effect, which is induced by the temperature of 
the gas is negligible in this case. 

By comparison, the Vishniac effect (1987) induces dis- 
tortions due to the correlations between the density and 
the velocity distributions in a fully ionized medium. The 
Vishniac effect is a second order one and it dominates on 
small angular scales (a few arcminutes or a few tens of 
arcminutes) . 
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We investigate here a first order eff'ect due to the bulk 

velocity of the ionized portion of the line of sight 

Observing the primordial anisotropies of the CMB is 
the only way to constrain the various models of structure 
formation and give information about the origin of the 
fluctuations. In fact, the scale-invariant spectrum of mass 
density is preferred in many models. The standard picture 
of structure formation in an inflationary model naturally 
creates fluctuations with such a spectrum and a gaussian 
statistics. Cosmic defects, such as strings, textures..., are a 
different source of fluctuations. Although they might give 
a scale- invariant spectrum, they introduce a non-gaussian 
characteristic to the statistics of the fluctuations. 

It is therefore important to determine whether or not 
the statistic of the observed primordial anisotropies is 
gaussian. In this context, it is thus also necessary to in- 
vestigate mechanisms leading to deviations from gaussian 
statistics in the standard model. We find that the fiuctu- 
ations induced by the early quasars are non-gaussian; this 
effect introduces a new component to the distortions and 
makes it harder to draw conclusions about the statistical 
nature of the primordial fluctuations. 

Quasar formation theory is still in its infancy. Observa- 
tions show the existence of quasars at redshifts as high as 
z Ri 5, although in hierarchical galaxy formation theories, 
the galaxies form late. We note that a recent work (Katz et 
al. 1994) shows that quasars could still be formed at high 
redshifts even in theories where galaxy formation occurs 
late such as in the hierarchical theories. The present work 
has no strong dependence on quasar formation theory, the 
presence of any high luminosity ionizing point sources will 
have the same kind of effect on the CMB, as long as it can 
roughly account for the reionization of the universe. 

For the cosmological parameters, we used h = 
Fo/lOOkms-^Mpc"! = 0.5, fio = 1 (flat universe) and 
rifc/i^ = 0.0135 fraction of baryonic matter; according to 
the nucleosynthesis constraints we have 0.01 < Ofc/i^ < 
0.015 (Walker et al. 1991). Hereafter, we assume that there 
are no heavy elements other than those created in the Big 
Bang (H, D, He, Li, ...). Moreover, we neglect the con- 
tributions other than those induced by hydrogen and he- 
lium, and we assume that the quasars emit their radiation 
isotropically, producing spherical ionized bubbles. Possi- 
ble geometrical effects were neglected in the present work. 
Owing to the high redshift range of interest, the hypoth- 
esis of homogeneity for the IGM gas is valid. 

These approximations are justified in view of other un- 
certainties on the reionization process. The goal of this 
paper is to show that a non negligible effect is induced by 
this process and not to compute a detailed model. 

In section 2 we discuss the relevant properties of an 
isolated bubble around a quasar and we compute the dis- 
tortion of the CMB induced by such a bubble. In section 
3 we compute the statistics over the whole sky of the dis- 
tortions induced. In section 4 we present the results and 
discuss them in relation to recent CMB measurements. 



2. Isolated early quasars and their effect on the 
CMB 

The specific quasar model we use is a good illustration for 
effects induced by other possible sources of ionization like 
young ultraluminous galaxies for example. 

We define the period of isolated quasars as the redshift 
range where most of the baryonic matter is still neutral 
while early quasars appear. 

An isolated quasar emits strong radiation - in partic- 
ular in the UV range - that can ionize the neutral inter- 
galactic gas around it at great distance, thus creating an 
ionized sphere. 

2.1. Ionizing radiation production rate 

We first compute the size of the ionized sphere surround- 
ing the quasar. Knowing the quasar specific luminosity in 
the B band, we deduce the emission rate S of UV photons 
ionizing the hydrogen atoms of the IGM. 

The intrinsic spectrum of individual quasars in the ion- 
izing UV domain is still uncertain but a good approxima- 
tion is to use a power law extrapolation of the optical flux. 

Several observations have been made in different spec- 
tral ranges (optical. X-rays ...) and constraint the values 
of the average spectral indexes in each range. A review 
of these observations was reported for example by Bech- 
told et al. (1987) who considered three spectral forms. We 
adopt the model consistent with the "medium" spectrum 
(Bechtold et al. 1987) which gives for the continuum spec- 
tral energy distribution radiated by quasars; 

La cx ly-^-'' (A < 1216 A). (1) 

The luminosity function for quasars defined in Sect. 3.1 
is written as a function of La the luminosity per spectral 
interval at the Lyman a frequency. It is therefore necessary 
to write S the UV emission rate - expressed in number of 
ionizing photons per unit time - as a function of La 

S = ^dv, (2) 

where hpi is the Planck constant and Uion is the ionizing 
frequency for hydrogen. 

Since the integral is dominated by the lower limit in 
frequency (ionizing limit for the hydrogen), the ionizing 
flux is not critically dependent on the assumed spectral 
shape. We get S = ALa where ^ is a constant de- 
pending on the normalization of the spectrum. We intro- 
duce the following notation: a characteristic luminosity 
^ai^c) = L*c which is independent on the quasar age 
and will be introduced in Sect. 3.1, and defining the ratio 
Iq = La/ LI, one gets 

S=^(^YlQ. (3) 

a = 1.5 is the spectral index of the quasar spectrum in 
the range of interest. 
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2.2. Temperature of the bubble 

If the intergalactic matter around the quasar - essentially 
composed of hydrogen and helium is ionized by the emit- 
ted radiation, the ionization of a hydrogen atom produces 
a photoelectron with an average kinetic energy defined by 
Ec = Ejjv — 13.6 eV, where Euv is the average energy of 
the ionizing photon. 

Since electrons and protons rapidly thermalize owing 
to the short time scale of the Coulomb interactions, the 
electrons transfer about half of their energy to the pro- 
tons, and the temperature of the plasma is thus given by 
T = Ec/Skb as long as the cooling by recombination is 
negligible; here kb is the Boltzmann constant. 

We can compute the average energy of the ionizing 
photon Euv = hpiV knowing the quasar spectrum and 
assuming that every photon is absorbed sooner or later: 2 10"^'* cm"^ s" 



(5) 



where rig is the electronic density rig = no(l + z)^ 
(no = 1.13 10~^ fifeft.^ cm~^ is the present electronic den- 
sity). This approximation is valid as long as the photoelec- 
tron energy does not exceed the photoionization energy 
that is to say that a photoelectron is unable to produce 
additional ionization. The recombination coefficient to 
all hydrogen levels is given by (Hummer & Seaton 1963) 



ar = 1.627 10"^^ T^"^/^ X 
(1 - 1.657 logT4 + 0.584 T^/^) cm-^s-\ 



(6) 



where T4 is the temperature of the ionized IGM in- 
side the Stromgren sphere in units of 10* K; thus: ~ 



(4) 



We find Euv ~ 40 eV, which leads to T « lO^K for the 
plasma temperature. A photoelectron even if it has an 
energy larger than 13.6 eV will have a very low probabil- 
ity to make a second ionization, the propagation of the 
ionization front being much faster than the electron ve- 
locity. Thus the number of ionization is very nearly equal 
to the number of ionizing photons. In that context, tak- 
ing into account the presence of helium is equivalent to 
a 10% increase of the baryon density with hydrogen only. 
Therefore, in the following, we use flbh^ = 0.015. 

In the bubbles surrounding the quasars, the gas is fully 
ionized, therefore there is no cooling due to neutral atoms. 
The ionized gas will cool by bremsstrahlung (free-free ra- 
diation), but, the cooling time is longer than the age of 
the universe. 

Likewise, we investigate the effect of the inverse Comp- 

ton cooling and find that it will become important for red- 
shifts greater than 10 h"^^^ (cooling time equal to the age 
of the universe), or for h = 0.5, z > 7.6. 

We check that the temperature of the bubble decreases 
by a factor 2 at 2; = 10, the Compton cooling does not 
change drastically the temperature. Since the exact tem- 
perature is not critical, we assume that it remains constant 
at T « 10^ K but undergoes the adiabatic cooling due to 
expansion. 

This approximation breaks down for redshifts signifi- 
cantly larger than 10. 

2.3. Radius and angular size 

We now compute the radius of the ionized region induced 
by the quasar and first the Stromgren sphere (stationary 
ionization bounded Hii region). 

The proper radius Eg of the Stromgren sphere is de- 
rived from 



Finally, the Stromgren radius is given by: 



Rs — 



35 



1/3 



(7) 



On the other hand, let us assume that a quasar turns on 



at a redshift 



and turns off at 



'off 



in the assumption 



of a lifetime tq short compared to the expansion time, we 
have 1 -|- Zon ~ 1 + Zof /. The proper radius R of an ionized 
bubble that is produced around it at a given time t (when 
it has not reached the stationary Stromgren sphere state) 
is computed by writing that every UV photon emitted by 
the quasar is absorbed through the photoionization of a 
neutral H atom; therefore the main quantity that rules the 
size of the bubble is the total number of ionizing photons 
that are emitted by a quasar during its lifetime: this is 
given by the product Stq. 

The proper radius R will reach a maximum value at 
the redshift z^ff (turn off of the quasar) Rmax{La, Zoff)', 
and then it undergoes the average Hubble expansion (the 
dynamical expansion due to the internal pressure of the 
bubble can be neglected). 

This maximum radius is the solution of: 



47r 



Stq. 



(8) 



At a given redshift z smaller than Zo//, the expansion 
leads to the radius given by the following relation: 



(9) 



defining R,, = {3AL*tq/4:TTnoY/'^ a characteristic radius 
associated with the characteristic luminosity, L* and the 
quasar lifetime tq we finally have: 



RilQ,z) = Rjlj\l + z)-\ 



(10) 



For tq ~ lO^yrs, we find for our choice of cosmological 



parameters Rc = 10 



\^0.015 J 



Mpc. 
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The comparison of the two computed radii (Eq. (7) and 

(10)) shows that for rcdshifts less than Zeq where 1 + Zeq = 
{tqaruoy^/^, (« 47 for VLbK^ = 0.0135 and tq = lO^rs) 
the ionized bubble is smaller than its Stromgren radius 
and therefore recombinations can be neglected. 

The angular size (radius) of the ionized bubble is 
6 = R/Da, where is the angular distance; for a flat 
universe it is given by 

Da = ^{l + z)-^[l + z-VT+^]. (11) 
■Ho 

Finally, by defining 9c = RcHq/2c, which is given by 9c ~ 

^{^y^'\^y^'\-h) arcmin, and using Eq. (10) 
we get 

9 = 9cl]i\l-{l + z)-'/^]-\ (12) 

which is proportional to {VLi,h^)~^/'^ h, in a flat imivcrsc. 
Figure 1 shows the angular radius of an ionized bubble 
as a function of redshift for different values of the quasar 
luminosity and lifetime. 



or = 6.65 10 cm^ is the Thomson cross section. We 
have thus: 

Tm = 2aTne{z)R{La,z), (14) 

putting Tc = <TTno{z)Rc the characteristic optical depth 
obtained for L*, we have 

Tm = Tclli\l + Zf. (15) 

The average optical depth over the spherical structure is 
Ta = (jTne{z)R{La, z) with R{La, z) = ^R{La, z), there- 
fore 

2 

Ta = -^Tm, (16) 

and it scales like {nbh^f/^ Q^^^^. Fi gure 2 shows the de- 
pendence of the bubble optical depth with both the life- 
time of the quasar and its luminosity, as a function of the 
redshift. One should notice in Fig. 2 that the typical op- 
tical depth of 10~^ (comparable to those of rich clusters) 
are reached for moderate redshifts. 




Fig. 1. Angular radius of an ionized bubble as a function of the 
redshift for three values of the quasar luminosity, where L* is a 
characteristic luminosity defined Sect. 3.1. The angular radius 
varies as L^^^ see Eq. (12). The series of thin Unes stands for 
the case of the angular radius computed with a quasar lifetime 
tq = 10^ yrs, whereas, the thick lines are for tq = 5 10^ yrs. The 
angular radius increases with (Eq. (12)). 



2.4. Optical depth 

Knowing the proper radius of the ionized sphere, we can 
compute the Thomson optical depth for the line of sight 
going through the center of the bubble: 

Tm = 2 axriedl, (13) 
Jo 



Fig. 2. Average optical depth of an ionized bubble versus 
redshift for different values of the luminosity with lifetime 
tq = yrs in the case of the thin lines. The thick lines 
represent the increase due to the increase of the lifetime 
tq = h 10' yrs. 



2.5. Thermal Sunyaev-Zeldovich effect 

Compton scattering of the CMB radiation on hot free elec- 
trons imprints on the incident photon spectrum a signa- 
ture known as the Sunyaev-Zeldovich thermal effect (SZ) 
as proposed by Zeldovich & Sunyaev (1969) and Sunyaev 
& Zeldovich (1972; 1980) by redistributing the photons in 
the spectrum with a general shift toward higher frequen- 
cies. This is an effect that has initially been proposed for 
clusters of galaxies, and it has already been observed in 
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the direction of several rich clusters (see for ex Birkinshaw 
1990; Birkinshaw ct al. 1991; Wilbanks ct al. 1994). 

The comptonization parameter, y, which characterizes 
the SZ distortion, depends on the bubble's characteristics 
and is essentially proportional to the line integral of the 
electron pressure along the line of sight. For an isothermal 
sphere, it is given by: 



y = 



(17) 



Te is the electronic temperature and rUeC^ is the electron 
rest mass. In the Raylcigh Jeans domain the distortion 
amounts to = —2y. It will be shown to be negligible 
because of the rather low temperature of the plasma in 
the bubble. 



2.6. Kinetic Sunyaev-Zeldovich effect 

Another type of interaction between the bubble and the 
CMB photons is the SZ kinetic effect (first-order Doppler 
effect). 

The natural reference frame to describe the SZ effect 
is the local standard of rest defined by the very large scale 
distribution of matter, i.c the frame in which the CMB 
is isotropic. If the bubble is at rest in this frame, then 
the thermal SZ effect (second order effect with respect to 
the average velocity of the electrons) is the only signature 
imprinted on the spectrum. 

The significance of the so called kinetic effect has been 
given by Sunyaev & Zeldovich (1972); the motion of a 
cluster with respect to the background leads (due to the 
Doppler effect) to an additional change of the radiation 
temperature in its direction, because of the finite optical 
depth associated to the bubble. The spectrum of the dis- 
tortion is indistinguishable from a primordial anisotropy 
spectrum which is equivalent to a simple temperature 
change. The observed radiation temperature is changed 
by St = AT/T = —{vr/c)T, in a direction crossing the 
bubble where r is the Thomson optical depth and is 
the radial component of the bubble peculiar velocity (a 
positive one corresponds to a recession velocity). 

In order to evaluate the relative temperature fluctu- 
ations due to the SZ kinetic effect associated with the 
spheres of plasma around early quasars or left by them 
after they stop emitting ionizing radiation, we need the 
optical depth (see Sect. 2.4) and the radial velocity of the 
bubbles. The average anisotropy associated with an ion- 
ized bubble is given by: 



(18) 



this gives us a direct relation between the luminosity 
of a quasar, its radial velocity, the redshift and 5t- 
As regards the cosmological parameters, we have 5t oc 

(06/l2)2/3j^-V6_ 



2. 7. Velocities 

In standard cosmological models, large scale velocities are 

assiimcd to be induced by gravity from mass density fluc- 
tuations: infall into over dense regions and outfall out of 
underdense regions. According to the linear gravitational 
instability theory, the rotational part of the velocity field 
is diluted by the expansion so that the present large scale 
cosmic velocity field is expected to be a potential flow. The 
velocity dispersion of the large scale structures is com- 
puted in the assumption of the linear theory; it is given 
by ay{z) = f7o„(l + z)-'^!'^ (Peebles 1980; 1993), where 

TOi; = (J^{z = 0). 

In the following, we will assume that the bubbles pro- 
duced by quasars are random substructures in the linear 
potential of large scale structures and thus in motion with 
them. The distribution of their radial velocities is therefore 
given by the large scale velocity field distribution. The ve- 
locities are dominated by the maximum of (Jot,(L), where 
L is the corresponding scale, the maximum takes place for 
scales about 10 Mpc which is still linear at redshifts of in- 
terest (z > 5). The velocity distribution is directly related 
to the distribution of primordial density fluctuations dur- 
ing the linear stage; in the assumption of gaussian statistic 
for the fluctuations, the peculiar one dimensional velocity 
distribution of the quasars and thus the ionized bubbles 
is also gaussian, with a velocity dispersion ay{z). 

Numerical simulations have been done to evaluate (Tod 
and showed for different cosmological models that it was 
in a range about 700km/s to 200km/s depending on the 
model parameters (Rhee 1993; Croft & Efstathiou 1994). 

The observed bulk flow within a sphere out to 6000 
km/s is 350 to 400 km/s (Faber et al. 1993). We take here 
a semi empirical approach using the z dependence given 
by the linear theory and a conservative value of (Tqi, = 
300 km/s. The coherence length of the bulk flow could be 
significantly larger than in the Cold Dark Matter (CDM) 
model and thus lead to a larger effect. 

The ratio of SZ thermal effect to kinetic effect is inde- 
pendent on the optical depth and given by: 



m,,,c^ Vr 



= 1.710" 



10^ K 



300 km/s 



(19) 



Hence the thermal effect is negligible compared with the 
kinetic effect. This is opposite to what happens for clus- 
ters, because the intracluster gas is much hotter (w 10^ K). 

Figure 3 shows the temperatiu'c fluctuation taking for 
the value of the velocity, the velocity dispersion, as a func- 
tion of redshift for different luminosities and lifetimes. 

3. Quasars distribution 

3.1. Quasars luminosity function 

The statistics of the measurable bubbles can be deduced 
if one knows the distribution and evolution of the ion- 
izing sources. This is done when we take into account 
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Fig. 3. Temperature fluctuation associated with an ionized 
bubble versus tlie redshift for a radial velocity given by the 
velocity dispersion (Sect. 2.7) and different quasar luminosities. 
Thin lines are obtained with quasar lifetime tq = lO' yrs and 
thick ones with tq = 5 10^ yrs. 



the quasars distribution given by their luminosity func- 
tion ^{La,z) which expresses for a given redshift z the 
comoving space density of quasars as a function of the 
luminosity. 

The quasars luminosity function is not well determined 
yet, in particular at high redshifts {z > 3), owing to the 
small number of observed sources. Therefore, we will con- 
sider different expressions for the luminosity evolution de- 
pending on the redshifts range of interest. 

According to recent surveys (Boyle, Shanks, & Peter- 
son 1988; Boyle, Jones & Shanks 1991), a pure luminosity 
evolution model - in which the comoving number density 
of quasars is assumed to be statistically constant with re- 
spect to the time while the luminosity varies with z - 
describes well the properties of the quasars population 
up to 2 « 3. For z > 3, following Warren et al. (1991) 
and Meiksin & Madau (1993) we introduce a pure den- 
sity evolution model to describe the population beyond 
z ^ 3 which states that the comoving number density de- 
cays exponentially by a factor of 2 per unit redshift (see 
Hartwick & Shade 1990 for discussion of the quasars lu- 
minosity function). 

The luminosity function can be written as: 



$(La,-2) = (l>{La,z)f{z,z) 



where 

(j){La,z) 



+ 



La 



(20) 



(21) 



is the best fit to the quasars luminosity function up to 
z « 2.9, (Boyle 1991), and / is a luminosity evolution 
function (see below). 



In a cosmological model with f2o = 1 and Hq = 

SOkms"-'^ Mpc"'^, the "best-fit" parameters for the model 
are (Boyle 1991) = 6.5 IQ-'^Mpc-^, /3i = 3.9, P2 = 1.5. 

The absolute B magnitude, Mb, variation with red- 
shift is expressed by Mb{z) = M^{z = 0) - 2.5kL log(l + 
z), where the B magnitude at the present epoch is 
Af^(z = 0) = -22.4 and fci, = 3.45 ("best-fit" param- 
eter). This defines a characteristic luminosity L^iz) = 

10-0.4Mb(.)+20.29 g^gg-l 

It finally gives 



La{0){l + z)''^ for z<Zc 
i* {zc) = L* otherwise. 



(22) 



where Zc is a characteristic redshift, Zc = 1-9 ("best- 
fit" parameter). The luminosity evolution models are ex- 
pressed by the function f{z,z) depending on the value of 
characteristic redshift z. In the constant comoving range 
for z < z = 3we have f{z, z) = 1, while in the exponential 
decay mode, for 2; > z = 3, f{z,z) = exp{—ix{z — z)); ^ 
is the decay constant. The factor 2 decline on the number 
of the observed quasars a high redshifts gives = 0.69. 

The comoving spatial density of the quasars in function 
of the redshift in the range of interest is therefore given 
by the following expression: 



(j){La,z) 



iop{-n{z - 3)) X 





La 


01 


La 




1 




+ 







if z> z = 3. 



(23) 



3.2. Quasars formation rate 

The luminosity function gives the number of the ob- 
served quasars but it can be associated with very different 
quasars production rate depending on their lifetimes. This 
is important because we need to include ionized bubbles 
created by already dead quasars. The number of bubbles 
at a given time is thus larger than the number of radiating 
sources at the same moment. Only under the assiimption 
that quasars have a lifetime much longer than the age of 
the universe, does the luminosity function represent the 
number of ionized spheres. The size of the bubble will in- 
crease faster than {l+z)~^ if the quasar lifetime is infinite. 

A more likely assumption is that the quasar lifetime 
is short relative to the age of the universe. In this case, 
although the quasar does not emit any more ionizing pho- 
tons, it has already produced a long lived bubble that 
undergoes expansion as we have seen in Sect. 2.3. 

We define the rate of production of quasars Q{La,z) 

as 



Q{La,z) 



HLa,z) 



t„ 



(24) 



where we assume that the luminosity function is constant 
during the short quasar lifetime tq. This quantity will thus 
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express the spatial density of ionized bubbles depending 
on the rcdshift and the luminosity of the quasars that have 
formed them. 

The theory of quasar formation and evolution is still 
very uncertain, therefore the quasar lifetime tq was alter- 
natively given two values tq = 10''yrs and tq = lO^yrs 
(Padovani et al. 1990). We note that decreasing the life- 
time, leads to the production of smaller and more numer- 
ous bubbles. 

3.3. Porosity parameter 

The Gunn-Petcrson (1965) test for hydrogen shows that 
the universe is ionized at redshift greater than 5. Re- 
cent measurements of the Gunn-Peterson effect for he- 
lium based upon the Hubble Space Telescope FOG obser- 
vations by Jakobsen et al. (1994), indicate that the Hen 
reionization of the universe is partial until z Ki 3 (Madau 
& Mciksin 1994). We suppose that the photoionization 
of the intergalactic medium is due to the production of 
Hii regions associated with UV sources (quasars); these 
regions expand and must overlap by Zion ~ 5. 

In order to evaluate the redshift Zion of the reionization 
of the universe for hydrogen, one can compute the porosity 
parameter P. It is the product of the volume associated 
with an individual ionized bubble by the spatial number 
density of the bubbles 



r+oo 
J z 



dt 



dzoff 



dzoff 



— R {La,z)(j){La,Zoff)dLo 



(25) 



thus be computed for these bubbles. Given the luminos- 
ity function and the rate of production of qiiasars we are 
able to give statistical properties for the relevant quanti- 
ties; we compute the number of bubbles as a function of 
the rcdshift, their angular size and the GMB temperature 
fluctuation they induce, the surface covering factor... 

The induced GMB temperature fluctuations will be 
dominated by a rcdshift range in which individual ionized 
bubbles induce a significant Doppler effect, this is when 
the hydrogen is still mostly neutral that is for redshifts 

Z ^ Ziorf 



10.0000 F 



U.O I 'yJU 

0.0010 
0.0001 



1 



Fig. 4. The porosity parameter is plotted {solid curve) as a 
function of the redshift. The horizontal dashed line stands for 
the reionization. Porosity parameter is independent on the life- 
time of the quasars. 



Since every UV photon must be absorbed in the IGM, 
the volume filling factor of ionized hydrogen is equal to 
the porosity parameter P (and not to 1 — exp(— P), which 
would be the case if the physics would fix the radius of 
the bubble rather than the ionized volume). Thus the com- 
plete reionization of the IGM occurs for P reaching 1. Fig- 
ure 4, which is plotted using the set of parameters given in 
Sect. 1, shows that the assumptions from previous sections 
lead to Zion ~ 5.7. Furthermore, the choice of the cos- 
mological parameters influence the value of the porosity 
parameter {P oc {ilbh^)~^ Qq^^^). The reionization 
redshift is independent on the quasar lifetime as long as 
the quasars production rate and associated lifetime leads 
to the observed luminosity function. Our assumptions are 
therefore compatible with the observed reionization of the 
universe. 

The integration over the luminosity takes into account 
the total contribution of the quasars to the reionization. 

At redshifts where P <C 1 , the universe is only partially 
ionized and the produced ionized bubbles are isolated in 
a neutral universe. All the quantities deflned Sect. 2 (an- 
gular size, optical depth and temperature fluctuation) can 



3.4- Number of bubbles 

Starting from the quasars production rate, we first com- 
pute the number of ionized bubbles N as a. function of 
their associated temperature fiuctuations St = AT/T, 
their angular sizes 0, redshift z, per units of solid angle. 
Knowing the number density of bubbles per units of vol- 
ume, N, which is given by: 



N 



/ Q{L„,z)dt. 
Jo 



(26) 



We derive the number of bubbles dN = NdV in the el- 
ementary volume dV = a^r"^ drdfl with adr = cdt and a 

the scaling factor a = {l + z)~^, dft is the viewing solid an- 

z)"^/^) assuming an Euclidean 



gle; with r 
universe. 



3cto(l-(l 



d'^N 
dzdO. 



p+oo 

/ *(ia, ^o//)(l + Zoff)-^/^dZoff, 
J z 



(27) 
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here k{z) = (1 + z)-^/^[l + (1 + z)-^/^f. 

Using the quantities computed in Sect. 2 (radius, an- 
gular size and optical depth of the ionized bubbles around 
quasars) which are function of the redshift and the quasar 
luminosity, wc can compute the number of bubbles at red- 
shift z with angular size 9 and absolute value of the tem- 
perature fluctuation \6t\ in the universe, per unit of solid 
angle as follow: 



dil,dzdhi9 din \6t\ 



\HoJ tc 



f{z)9\dT\x 



Hh)9{vr, z) d[lQ - Iq{0, z)] dlQ 



5[5t - 6T{Vr,lQ,Z 



r+oo 

)]dVr / J 
J z 



p{zoff)dzoff; 



(28) 



where 6 is the Dirac function, g{vr, z) is the gaussian dis- 
tribution for the radial peculiar velocities of the quasars, 

/(^) = (l+z)-9/2[l + (l + z)-V2], 

Pilots) = (1 + ZoSfY'^ exp[-0.69(^o// - 3)] 
and Iq{6,z) is derived from Eq. (12) and Vr from Eq. 
(18). 

There are as many negative sources as positive ones, 
therefore the factor 2 is introduced when we take into 
account the negative and positive radial velocities of the 
bubbles and give the counts as a function of We note 
that in a flat universe, the number of bubbles scales like 

The results have been numerically computed and are 
plotted Fig. 5 and Fig. 6, integrated over the redshifts 
{z > 6) for different values of quasar lifetime {tq = 10* yrs 
and tq = lO^'yrs). 

Next, we also compute the covering surface parame- 
ter Cs which is the fraction of the sky covered by ionized 
bubbles inducing a distortion of amplitude 6t and angular 
size 6 at redshift z. This is done by multiplying the previ- 
ous counts by the solid angle occupied by the bubbles; we 
obtain: 



Cs = 2 X 6 



Hr 



to 



fiz) X 



J J TTe'^^{lQ)g{vr,z)dlQdvr5[lQ-lQ{e,z)] 

r+oo 

6[5t - 6T{Vr,lQ,z)] / l{Zoff)dZoff. 

J 2 



(29) 



Factor 2 takes into account the sources with both negative 
and positive velocities, and Cs oc (Oft/i^)"^/^ h~'^ for a flat 
universe. 

The results plotted Fig. 7 and Fig. 8 are also integrated 
over z and given for the two previous values of the quasar 
lifetime tq. 



dN/dln(0)/dln((5T/T) per sr 
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Integrated over z>6 with tq=le7 yrs 



Fig. 5. The number of ionized bubbles per unit solid angle, 
logarithmic interval of \6t\ logarithmic interval of the angular 
radius 6 is plotted as a function of the temperature fluctua- 
tion log^T and log^ in arcmin. This number is integrated over 
the redshifts {z > 6) and computed for lifetime tq = lO'^yrs 
according to Eq. (28). 



dN/dln(0)/dln((5T/T) per sr 
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-6.0 



-4.5 



-5.5 -5.0 
log(<5T/T) 

Integrated over z>6 with tq=le8 yrs 



-4.0 



Fig. 6. The same quantity is plotted for a quasar lifetime 
tq = 10^ yrs. The bubbles are greater but less numerous. 



3. 5. Anisotropics during fully ionized period 

When the IGM is totally ionized at ^; < Zion = 5.7, 
wc evaluate the rms temperature fluctuation. The optical 
depth is almost constant in this redshift range (z < Zion) 
whereas the velocities are variable. Therefore, we compute 
the average of the temperature fluctuation along a line of 
sight 



to 



HpVr dt. 



(30) 



The fluctuation is due to the peculiar velocity variation 
which is expressed by the velocity dispersion. In the lin- 
ear perturbation theory, it is shown that the velocities de- 
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Surface covering factor 
2.0 I 

S 1.5- 

6 



O 




-6.5 -6.0 -5.5 -5.0 -4.5 -4.0 
log((5T/T) 

Integrated over z>6 with tq=le7 yrs 

Fig. 7. The surface covering factor of the ionized bubbles per 
logarithmic interval of \5t\ and logarithmic interval of the an- 
gular radius 8 is plotted as a function of the temperature fluc- 
tuation (log^r) and logd in arcmin. It is integrated over the 
redshifts {z > 6) and computed for lifetime tq = lO'^yrs ac- 
cording to Eq. (29) 

dN/dln(0)/dln(5T/T) per sr 
2.0 



.S 1.5- 




-6.5 -6.0 -5.5 -5.0 -4.5 -4.0 
log(5T/T) 

Integrated over z>6 with tq=le8 yrs 

Fig. 8. The surface covering factor is plotted for quasar lifetime 

tg = 10* yrs 

pond upon the scales of the density perturbation that are 
associated with them and it is also shown that in a CDM 
model the coherence length is short. Other cosmological 
models could give a stronger effect. 

Using the velocity dispersion of the distribution, we 
compute the rms temperature fluctuation {ST)rms when 
the IGM is fully ionized, and find 

(«™»1.5I0-«(5^)(A)". (31) 

We note that the effect due to the velocity variation is 
one order of magnitude lower than the rms fluctuation 
detected by COBE (Mather et al. 1994) and it is at the 



limit sensitivity of the modern detectors, the average effect 

of velocity field in the fully ionized universe is therefore 
weaker than the effect for a partially ionized universe. 

4. Discussion 

We discuss here the effects of ionization of the IGM by 

early formed quasars on the cosmic microwave background 
as calculated in previous sections. 

We have shown the thermal SZ effect was negligible 
compared to the kinetic one. The amplitude of the largest 
temperature fiuctuations induced by the latter mechanism 
are comparable to the the primordial fluctuations at scales 
comparable with the scales of the Doppler peaks (10 to 
1°). These secondary fiuctuations have no spectral signa- 
ture to distinguish them from the primordial ones. Hence, 
future observations and detections of the CMB at small 
scales with sensitivity about 10~® will have to take into 
account the possible presence of such objects in the sky. 

When the porosity parameter is lower than 1 for z > 
Zion, we have computed the number of ionized bubbles per 
unit of solid angle, with respect to their angular radius and 
to St the temperature fluctuation they induce integrating 
over the redshifts. 

We plot (Fig. 6) the number of kinetic effect sources 
integrated over redshifts (z > 6) per unit solid angle, first 
for lifetime tq = 10* yrs. We find that there are approx- 
imately 10^ objects per unit solid angle with associated 
temperature fiuctuations some 10~^ and angular radius 3 
to 10 arcmin randomly distributed over the sky and about 
10 ionized bubbles per unit solid angle with temperature 
fluctuations with about 3 10~^ and radius about 12 to 30 
arcmin. With the same set of parameters, we plot (Fig. 8) 
the surface covering factor integrated over redshifts versus 
the amplitude of the anisotropics and their angular radius. 

The counts depend significantly on the quasar lifetime 
value, in fact, it was shown that decreasing the lifetime tq 
decreases the size of the bubbles, but on the other hand, 
it increases the rate of formation of quasars and thus gives 
more numerous bubbles. Figures 5 and 7 show respectively 
the number of bubbles per unit solid angle and covering 
surface versus St and the angular radius for tq = 10*" yrs. 
In this configuration, we find that there are 100 objects per 
steradian with St about 10~^ and radius 5 to 12 arcmin. 

Therefore, early quasars will induce distortions with 
different statistical properties than the primordial ones 
and produce rare fluctuations with amplitude of up to 
10"^. Although the rms value of the fluctuations induced 
by the inhomogeneous reionization is smaller, by an order 
of magnitude, than the rms of the fiuctuations detected 
by COBE (COsmic Background Explorer) (Smoot et al. 
1992; Mather et al. 1994), the largest ones at small angu- 
lar scales could be dominated by the effect discussed here. 
Hence high accuracy measurements of the power spectrum 
of the CMB anisotropies such as those considered for fu- 
ture space experiments (better than 1%) and search for 
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deviations to Gaussian statistics must take into account 

this kind of sources. We find that most of the anisotropics 
with high level AT/T and small angular scales are given 
by quasars at high redshifts. 

Observations from the first flight of the Medium Scale 
Anisotropy Measurement (MSAM) (Cheng et al. 1994) 
with a 30' beam instrument, showed the presence of 
two unresolved sources from a survey of 7.610^''sr with 
5t ~ 3 10~^ consistent spectrally with CMB anisotropies 
; this implies about 200 of such sources per steradian. 

The temperature fluctuations produced by the eflfect 
that we are studying - ionization by early quasars - are 
indistinguishable from primordial anisotropies. Further- 
more, the computations done for our model with short 
quasar lifetime tq = lO^yrs predicts about 100 events 
per unit solid angle with angular diameter about 20' and 
6t ~ 310"^, consistent with the results of the MSAM 
observations. In this context, the statistics of such events 
on the sky will be an important criterion to attribute an 
origin to them. 

Figure 7 - representing the surface covering factor - 
shows that the ionized bubbles around quasars are close 
to cover the whole sky although they do not fill the space 
until z < Zion = 5.7. Therefore, we can no longer con- 
sider the induced temperature fluctuations as individual 
sources. On the contrary, it is necessary to get statisti- 
cal informations about the distribution of the anisotropies 
such as the standard deviation, mean value, kurtosis... 
Hence we have simulated maps of the anisotropies induced 
by the quasars. 

The maps were drawn under simple geometrical as- 
sumptions (spherical bubbles), using the counts of ionized 
bubbles associated with quasars characterized by an angu- 
lar radius and a temperature fluctuation given in Eq. (28) 
and using the fact that there is a simple relation between 
the angular radius of the bubble and the luminosity of the 
quasar that induced it, Eq. (12). We put a limit on the 
number of bubbles; for numbers of bubbles greater than 
this limit, a gaussian drawing is done otherwise we make 
a poissonnian drawing. The positions of the bubbles over 
the map are random, this is a rather good assumption 
since the clustering of quasars is still uncertain. 

We check that the standard deviation, mean, maxi- 
mum and minimum values are found consistent with the 
bubbles counts. But as could be expected, the distribu- 
tion of the temperature fluctuations induced by quasars 
is far from being gaussian. Indeed, for a typical value of 
the quasar lifetime tg = lO^yrs we flnd that the kurto- 
sis is equal to 1.4. This is an important result that might 
be taken into account if one wants to test the gaussian- 
nity of the primordial anisotropies of the CMB, to prove 
whether or not the primordial fluctuations arc adiabatic 
(thus due to quantum fluctuations of the vacuum energy) 
for which the statistic is gaussian or to topological de- 
faults (strings, domain walls, ...). In fact, the presence of 



early ionizing sources (quasars here) induces additional 
non gaussian anisotropies with no spectral signature. 

Cosmological models show that there is an angular cut 
off in the spectrum of the primordial anisotropies. The 
bubbles contribution to the observed CMB at small an- 
gular scales could be important or may be dominant, in 
particular at scales where the primordial fluctuations are 
"washed out" because of the damping associated with the 
width of the the last scattering surface. This effect is sug- 
gested Fig. 9, where we plot both the power spectra of 
the primordial temperature fluctuations and the bubbles 
in terms of the spherical harmonics coefficients (see for 
example White, Scott & Silk 1994); taking tg = lO'^yrs. 
Our results show that the secondary anisotropies of the 




10 100 1000 

Spherical harmonics order (I) 



Fig. 9. Power of the anisotropies induced by the bubbles and 

the primordial fluctuations in the spherical harmonics versus 
the multipoles. Solid curve for the primary CMB and dashed 
curve for the fluctuations induced by the quasars we take 
tg = W'' yrs 

CMB, produced by the kinetic SZ effect in ionized bub- 
bles around quasars, must be taken into account. 

In the future, observations along the lines of sight 
of known quasars at high redshift showing the Gunn- 
Peterson effect, could be an interesting test to detect this 
kinetic effect induced by the ionized bubbles. The follow- 
ing calculations are possible only for the bubbles associ- 
ated with the strongest temperature fluctuations for which 
primordial anisotropies do not perturb too much the ob- 
servations. In that case, we could be able to derive the 
bubble peculiar velocity. 

In Sect. 2.6, We have derived a simple relation be- 
tween the temperature fluctuation associated with an ion- 
ized bubble, the radial peculiar velocity and the luminos- 
ity of the quasar that has produced it Eq. (18). On the 
other hand, Eq. (12) gives the the expression of the an- 
gular radius of the the bubble as a function of the quasar 
luminosity. Therefore, if we observe a bubble inducing a 
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temperature fluctuation, we directly derive the luminosity 

from the angular size and knowing this we have a direct 
measure of the peculiar radial velocity. The latter is given 
by the following expression 

where 6c and Tc are defined respectively in Sect. 2.3. and 
Sect. 2.4. 

In the hypothesis of an ionization due to individual 
sources of radiation, the whole formalism we develop re- 
mains valid whatever the kind of sources, early galaxies, 
first generation of stars, ... The efi'ect of the presence of 
sources of luminosity of the same order as the quasars on 
the CMB will be the same as the one described for quasars. 
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